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LONG-TERM  GOALS 

To  develop  “exaet”  numerieal  methods  and  visualization  teehniques  for  the  study  of  propagation  of 
aeoustie  energy  in  shallow  water  in  the  time  domain.  Exaet  in  this  eontexts  means  the  methods  plaee  no 
restrietions  on  the  underlying  physies  of  the  environment. 

OBJECTIVES 

To  develop  new,  and  enhanee  existing,  numerieal  methods;  to  establish  the  aeeuraey,  robustness,  flexibil¬ 
ity,  and  traetability  of  these  methods;  and  to  apply  these  methods  to  meaningful  praetieal  problems.  To 
ereate  a  robust  software  suite  for  the  applieation  of  the  methods  and  to  develop  teehniques  for  visualizing 
the  propagation  of  sealar  and  veetor  fields  in  eomplex  environments. 

APPROACH 

We  have  foeused  our  attention  on  the  development  and  use  of  finite-differenee  time-domain  (EDTD) 
methods.  These  time-domain  teehniques,  whieh  are  flexible,  robust,  and  generally  simple  to  implement, 
have  been  used  used  extensively  in  the  past  by  the  eleetromagneties  eommunity  to  solve  a  wide  range  of 
problems.  We  have  established  the  aeeuraey  of  our  proposed  aeoustie  EDTD  teehniques  by  eomparing 
numerie  results  to  exaet  solutions  for  eanonieal  problems  or  to  results  obtained  using  other  numerieal 
methods.  The  EDTD  method  is  eomputationally  intensive  and,  henee,  to  inerease  the  ease  with  whieh 
EDTD  ean  be  applied  to  large  problems,  we  are  developing  eode  whieh  ineorporates  the  Parallel  Virtual 
Maehine  paradigm  to  permit  the  solution  of  a  problem  using  a  eluster  of  workstations. 

WORK  COMPLETED 

Eundamental  aspeets  of  material  boundaries  in  EDTD  models  were  studied  and  quantified;  inherent  prop¬ 
erties  and  artifaets  of  the  EDTD  grid  were  rigorously  explained;  EDTD  algorithms  were  developed;  the 
eorresponding  eomputer  programs  were  written;  and  information  was  disseminated  via  journal  publiea- 
tions,  eonferenee  presentations,  and  the  Web. 
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RESULTS 


The  FDTD  method  is  obtained  by  diseretizing  the  differential  equations  that  govern  the  underlying  sys¬ 
tem.  Using  a  Cartesian  grid,  the  method  provides  an  exeeedingly  simple  way  in  whieh  to  express  future 
fields  (i.e.,  unknown  fields)  in  terms  of  past  fields  (known  fields).  For  propagation  in  a  homogeneous 
region,  the  traditional  FDTD  method  is  aeeurate  to  seeond-order — that  is,  doubling  the  number  of  grid 
points  per  wavelength  reduees  inherent  numerieal  errors  by  a  faetor  of  four.  However,  the  behavior  and 
aeeuraey  of  fields  at  material  interfaees  are  mueh  more  eomplieated. 

We  previously  derived  exaet  expressions  for  the  transmission  and  refleetion  eoeffieients  for  fields  nor¬ 
mally  and  obliquely  ineident  on  planar  boundaries  and  showed  that  the  behavior  of  the  fields  on  one  side 
of  the  boundary  depends  on  the  diseretization  on  the  other  side  [1,2].  However  this  work  did  not  address 
the  more  ehallenging  problem  of  the  errors  eaused  by  the  inherent  “stair-step”  approximation  of  material 
boundaries  that  are  eontinuously  varying.  To  date  there  has  been  no  quantifieation  of  these  errors  for  an 
arbitrary  boundary.  Nevertheless,  for  a  given  objeet  one  ean  eodify  these  errors  through  measurements. 
We  have  used  measurements  to  demonstrate  the  aeeuraey  of  the  FDTD  method  when  applied  to  various 
problems  ineluding  randomly  rough  pressure-release  surfaees  [3,4,5],  penetrable  surfaees  [6],  and  rigid 
boundaries  [7,8]. 

In  all  our  work  modeling  eontinuously  varying  boundaries  we  have  developed  and  employed  loeally- 
eonformal  sehemes  whieh  use  modified  versions  of  the  FDTD  updated  equations  for  nodes  adjaeent 
to  material  interfaees  but  whieh  use  the  standard  equations  elsewhere.  In  [5,7,8]  we  showed  the  im¬ 
provements  that  sueh  loeally-eonformal  sehemes  ean  yield  while  introdueing  essentially  no  additional 
eomputational  eost.  To  demonstrate  the  improvement  that  ean  be  realized  using  a  loeally-eonformal 
seheme,  eonsider  the  rigid  sphere  shown  in  Fig.  1.  The  figure  shows  the  diseretized  version  of  the  sphere 
as  it  would  appear  in  a  Cartesian  grid  when  diseretized  using  a  radius  of  eight  eells.  The  sphere  is  in- 
sonified  with  a  broad-speetrum  pulsed  plane  wave  and  the  seattered  field  is  measured  along  the  heavy 
line  shown  in  Fig.  1.  The  time-domain  results  ean  be  Fourier  transformed  to  obtain  the  seattered  field  at 
any  frequeney  for  whieh  the  ineident  pulse  has  suffieient  speetral  eontent.  The  results  in  Fig.  2  show  the 
seattered  field  eorresponding  to  a  wavelength  whieh  is  modeled  using  9.853  points  per  wavelength.  The 
exaet  result  is  plotted  together  with  those  obtained  using  the  FDTD  seheme  with  a  uniform,  or  stairstep, 
grid  (UFDTD)  and  using  a  loeally-eonformal  seheme  (CFDTD).  Clearly  the  CFDTD  approaeh  provides 
a  signifieant  improvement  over  UFDTD.  Ten  points  per  wavelength  is  often  the  eoarsest  diseretization 
employed  in  FDTD  simulations.  For  the  simulation  geometry  shown  here  the  agreement  between  simu¬ 
lations  and  the  exaet  result  improves  as  the  diseretization  is  inereased,  but  regardless  of  the  diseretization 
the  loeally-eonformal  CFDTD  method  provides  results  superior  to  those  of  UFDTD. 

Despite  the  faet  that  we  have  used  a  loeally-eonformal  seheme  in  the  past  to  model  penetrable  objeets 
[6],  the  improvement  realized  by  using  the  loeally  eonformal  seheme  was  typieally  small.  We  have  re- 
eently  eompleted  a  rigorous  and  systematie  eomparison  of  various  low-eost,  loeally-eonformal  sehemes 
for  penetrable  objeets.  Various  sound  speeds,  densities,  objeet  sizes,  and  diseretizations  were  used.  Un¬ 
fortunately  there  was  no  elear  “winner”  in  terms  of  the  best  seheme.  However,  when  eonsidered  over 
all  the  test  eases,  the  seheme  whieh  typieally  performed  the  best  was  arguably  the  simplest  and  elosely 
akin  to  the  staireasing  teehnique:  densities  should  be  assoeiated  with  pressure  nodes,  and  veloeity  nodes 
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Figure  1:  Representation  of  the  3D  computational  domain  used  to  obtain  the  scattered  pressure  from  a 
rigid  sphere.  The  sphere  is  depicted  in  accordance  with  its  stair-step  representation.  The  incident  plane 
wave  travels  in  the  +z  direction  and  the  scattered  pressure  is  collected  along  the  position  of  the  heavy 
vertical  line. 


Figure  2:  Scattered  pressure  at  9.853  points  per  wavelength.  The  independent  variable  is  position  in 
terms  of  cells  where  zero  corresponds  to  the  same  z  location  as  the  center  of  the  sphere. 
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whose  neighboring  pressure  nodes  span  the  material  boundary  should  employ  the  average  density  in  their 
update  equations.  This  work  will  be  submitted  soon  to  an  appropriate  journal. 

We  continue  to  explore  several  new  implementations  of  the  FDTD  method  (proposed  by  others)  which 
seek  to  minimize  dispersive  and  anisotropic  errors  inherent  in  all  2-  and  3-D  FDTD  schemes.  We  have 
developed  comparisons  that  provide  insight  into  the  techniques  that  are  not  easily  garnered  from  the  pub¬ 
lications  in  which  they  were  originally  presented.  Some  of  this  work  will  appear  in  IEEE  Transactions 
on  Microwave  Theory  and  Techniques  [9]  and  was  presented  at  the  2002  URSFAntennas  and  Propaga¬ 
tion  Symposium  [10].  Notably,  we  demonstrated  that  many  of  the  wavelet-based  schemes,  which  have 
attracted  some  advocates,  are  not  superior  to  an  FDTD  scheme  that  uses  the  same  spatial  stencil  and  the 
same  “computational  effort”  (i.e.,  operations  per  a  given  temporal  advancement  of  the  fields). 

Along  these  same  lines,  we  have  adapted  a  variation  of  the  promising  FDTD  scheme  proposed  by  Eric 
Forgy  {IEEE  Transactions  on  Antennas  and  Propagation,  50(7):983-996,  2002).  This  algorithm  suffers 
much  less  grid  dispersion  and  anisotropy  than  more  traditional  FDTD  formulations  but  still  retains  the 
local  nature  of  the  standard  update  equations.  The  acoustic  implemention  will  be  described  in  a  future 
publication. 

It  has  long  been  known  that  the  Yee  FDTD  algorithm  can  provide  exact  solutions  to  one-dimensional 
problems  when  operated  at  the  so-called  magic  time  step.  Here  “exact”  is  taken  to  mean  the  field  prop¬ 
agates  without  dispersion  error  or  other  numeric  artifacts.  Unfortunately  there  is  no  magic  time  step  in 
higher  dimensions.  However  we  have  recently  developed  a  theoretical  framework  for  multi-dimensional 
algorithms  that  have  the  same  exact  properties  as  the  one-dimensional  Yee  algorithm  when  operated  at 
the  magic  time  step.  The  proposed  technique  uses  vector  operators  which,  instead  of  being  defined  at  a 
point  such  as  with  the  usual  gradient,  divergence,  and  curl  operators,  are  defined  over  spheres.  Due  to 
their  inherent  symmetry,  these  spatial  operators  have  the  same  properties  in  all  directions.  With  a  judi¬ 
cious  choice  of  the  temporal  step  size  the  temporal  errors  can  cancel  the  spatial  errors  and  the  algorithm 
is  exact.  However,  although  the  framework  for  the  algorithm  has  been  developed,  no  practical  (i.e.,  com¬ 
putationally  efficient)  algorithm  has  yet  been  developed.  Nevertheless  proof-of-concept  implementations 
of  the  algorithm  (which  are  quite  computationally  expensive)  have  been  used  to  demonstrate  the  validity 
of  the  technique  and  the  improvements  the  algorithm  can  provide  over  other  FDTD  implementations. 
The  algorithm  also  has  interesting  properties  such  as  unconditional  stability  for  an  arbitrary  temporal 
step  size.  We  are  continuing  to  explore  this  algorithm  and  will  be  publishing  our  findings  in  the  near 
future. 

We  continue  to  maintain  a  Web  site,  www .  f  dtd .  org,  that  seeks  to  list  all  archival  publications  related 
to  the  FDTD  method.  This  site  solicits  input,  in  the  form  of  comments  posted  about  work  appearing  in 
the  archival  literature,  from  the  entire  community  interested  in  the  FDTD  method  (whether  applied  to 
acoustics,  electromagnetics,  or  solid  mechanics). 

Our  investigations  of  the  discretized  worlds  of  FDTD  methods  have  led  us  to  a  better  understanding 
of  numeric  artifacts  associated  with  resonances  and  to  ways  of  alleviating  these  artifacts.  Part  of  this 
work  was  presented  as  an  invited  talk  in  a  special  session  organized  by  Prof.  Allen  Taflove  (one  of  the 
co-founders  of  the  FDTD  method)  [11].  This  work  is  further  described  in  a  paper  submitted  to  IEEE 


4 


Transactions  on  Antennas  and  Propagation  [8]. 

A  recent  publication  by  John  Pendry  (Phys.  Rev.  Lett,  85:3966,  2000)  which  described  the  use  of 
backward-wave  (BW)  materials  to  make  a  “perfect  lens”  caught  our  attention.  BW  materials  are  dis¬ 
persive  materials  whose  direction  of  phase  propagation  is  antiparallel  to  the  direction  of  power  flow.  BW 
behavior  can  exist  in  both  acoustic  and  electromagnetic  systems.  We  were  interested  in  modeling  BW 
material  using  the  FDTD  method.  Our  initial  results  were  not  consistent  with  those  one  would  expect 
from  an  initial  inspection  of  the  theory.  Eventually  we  discovered  that  the  dual,  offset  grids  employed 
in  the  FDTD  method  (i.e.,  the  dual  pressure  and  velocity  grids  or  the  dual  electric  and  magnetic  field 
grids)  can  introduce  significant  numeric  artifacts  when  modeling  BW  materials.  The  offset  in  the  grids 
can  introduce  a  boundary  layer  that  has  the  material  properties  of  neither  the  BW  material  nor  the  sur¬ 
rounding  medium.  Our  investigations  were  presented  in  an  invited  talk  at  the  2002  URSEAntennas  and 
Propagation  Symposium  [13]  and  in  a  paper  which  was  published  in  Physical  Review  Letters  B  [14]. 

IMPACT/APPLICATIONS 

Accurate  and  flexible  numerical  methods  allow  researchers  to  conduct  any  number  of  experiments  with¬ 
out  having  to  resort  to  actual  field  experiments,  i.e.,  the  experiment  is  conducted  on  the  computer.  Al¬ 
though  numerical  methods  will  never  supplant  field  experiments,  numerical  methods  (when  used  within 
their  “region  of  validity”)  do  provide  an  extremely  cost-effective  means  of  conducting  controlled  ex¬ 
periments.  Our  work  enables  more  accurate  and  more  efficient  numerical  solutions  to  a  wide  range  of 
problems  in  acoustics,  electromagnetics,  and  continuum  mechanics. 

TRANSITIONS 

Much  of  the  knowledge  we  have  gained  has  been  disseminated  via  publications  and  conference  presen¬ 
tations.  Additional  material  is  available  via  the  Web.  Please  refer  to  the  Web  site  given  in  the  header  for 
copies  of  the  Pi’s  publications  or  www .  f  dtd .  org  for  other  material  pertinent  to  the  FDTD  method. 

RELATED  PROJECTS 

This  work  is  related  to  research  in  both  high-frequency  acoustics  and  long-range  propagation.  Numerical 
models,  such  as  the  FDTD  method,  can  be  used  to  predict  the  fields  scattered  from  small  objects  under 
short- wavelength  insonification  or  the  propagation  of  long-wavelength  signals  over  limited  regions  of 
the  ocean.  Additionally,  this  work  is  related  to  several  other  ONR-sponsored  researchers  including  Shira 
Broschat,  Eric  Thorsos,  Philip  Marston,  and  Ralph  Stephen. 

REFERENCES 

[1]  J.  B.  Schneider  and  R.  J.  Kruhlak,  “Plane  Waves  and  Planar  Boundaries  in  FDTD  Simulations,”  IEEE 
AP-S  International  Symposium  and  URSI  Radio  Science  Meeting,  Salt  Fake  City,  UT,  Jul.  2000. 

[2]  J.  B.  Schneider  and  R.  J.  Kruhlak,  “Dispersion  of  Homogeneous  and  Inhomogeneous  Waves  in  the 
Yee  Finite-Difference  Time-Domain  Grid,”  IEEE  Trans.  Microwave  Theory  and  Techniques,  vol.  49,  no. 
2,  pp.  280-287,  2001. 


5 


[3]  F.  D.  Hastings,  J.  B.  Schneider,  and  S.  L.  Broschat,  “A  Monte  Carlo  FDTD  Technique  for  Rough 
Surface  Scattering,”  IEEE  Trans.  Antennas  Propagat.,  vol.  43,  no.  11,  pp.  1183-1191,  1995. 

[4]  F.  D.  Hastings,  J.  B.  Schneider,  and  S.  L.  Broschat,  “A  Finite-Difference  Time-Domain  Solution  to 
Scattering  from  a  Rough  Pressure-Release  Surface,”  J.  Acoust.  Soc.  Am.,  vol.  102,  no.  6,  pp.  3394-3400, 
1997. 

[5]  J.  B.  Schneider,  C.  L.  Wagner,  and  R.  J.  Kruhlak,  “Simple  Conformal  Methods  for  FDTD  Modeling 
of  Pressure-Release  Surfaces,”  J.  Acoust.  Soc.  Am.,  vol.  104,  no.  6,  pp.  3219-3226,  1998. 

[6]  F.  D.  Hastings,  J.  B.  Schneider,  S.  L.  Broschat,  and  E.  I.  Thorsos,  “An  FDTD  Method  for  Analysis 
of  Scattering  from  Rough  Fluid-Fluid  Interfaces,”  IEEE  Journal  of  Oceanic  Engineering,  vol.  26,  no.  1, 
pp.  94-101,2001. 

[7]  J.  B.  Schneider  and  J.  G.  Tolan,  “A  Locally  Conformal  Method  for  Modelling  Rigid  Boundaries  in 
the  FDTD  Method,”  J.  Acoust.  Soc.  Am.,  vol.  108,  no.  5,  pt.  2,  pp.  2563,  Newport  Beach,  CA,  Dec.  2000. 

[8]  J.  G.  Tolan  and  J.  B.  Schneider,  “Locally  Conformal  Method  for  Acoustic  Linite-Difference  Time- 
Domain  Modeling  of  Rigid  Surfaces,”  submitted  to  J.  Acoust.  Soc.  Am. 

[9]  K.  L.  Shlager  and  J.  B.  Schneider,  “Comparison  of  the  Dispersion  Properties  of  Several  Low- 
Dispersion  Linite-Difference  Time-Domain  Algorithms,”  to  appear  in  IEEE  Trans.  Antennas  Propagat., 
Dec.,  2002. 

[10]  K.  J.  Shlager  and  J.  B.  Schneider,  “A  2-D  Dispersion  Analysis  of  the  W-MRTD  Method  Using  CDL 
Biorthogonal  Wavelets,”  IEEE  AP-S  International  Symposium  and  URSI  Radio  Science  Meeting,  vol.  3, 
pp.  244-247,  San  Antonio,  TX,  Jun.  2002. 

[11]  C.  L.  Wagner  and  J.  B.  Schneider,  “Using  the  Dispersion  Relation  to  Understand  Linite-Difference 
Time-Domain  Worlds,”  International  Conference  on  Electromagnetics  in  Advanced  Applications  (ICEAA 
01),  Torino,  Italy,  pp.  375-378,  Sep.  2001. 

[12]  C.  L.  Wagner  and  J.  B.  Schneider,  “On  the  Analysis  of  Resonators  Using  Linite-Difference  Time- 
Domain  Techniques,”  submitted  to  IEEE  Trans.  Antennas  Propagat. 

[13]  M.  W.  Leise,  P.  J.  Bevelacqua,  and  J.  B.  Schneider,  “Backward- Wave  Meta-Materials  for  Perfect 
Lenses,”  IEEE  AP-S  International  Symposium  and  URSI  Radio  Science  Meeting,  San  Antonio,  TX,  Jun. 
2002. 

[14]  M.  W.  Eeise,  P.  J.  Bevelacqua,  and  J.  B.  Schneider,  “Effects  of  Surface  Waves  on  the  Behavior  of 
Perfect  Lenses,”  Physical  Review  B,  vol.  66,  no.  3,  035 113  (five  pages),  2002. 

PUBLICATIONS 

[1]  J.  G.  Tolan  and  J.  B.  Schneider,  “Locally  Conformal  Method  for  Acoustic  Linite-Difference  Time- 
Domain  Modeling  of  Rigid  Surfaces,”  submitted  to  J.  Acoust.  Soc.  Am. 


6 


[2]  K.  L.  Shlager  and  J.  B.  Schneider,  “Comparison  of  the  Dispersion  Properties  of  Several  Low- 
Dispersion  Finite-Difference  Time-Domain  Algorithms,”  to  appear  in  IEEE  Trans.  Antennas  Propagat., 
Dec.,  2002. 

[3]  K.  J.  Shlager  and  J.  B.  Schneider,  “A  2-D  Dispersion  Analysis  of  the  W-MRTD  Method  Using  CDF 
Biorthogonal  Wavelets,”  IEEE  AP-S  International  Symposium  and  URSI  Radio  Science  Meeting,  vol.  3, 
pp.  244-247,  San  Antonio,  TX,  Jun.  2002. 

[4]  C.  E.  Wagner  and  J.  B.  Schneider,  “On  the  Analysis  of  Resonators  Using  Einite-Difference  Time- 
Domain  Techniques,”  submitted  to  IEEE  Trans.  Antennas  Propagat. 

[5]  M.  W.  Eeise,  P.  J.  Bevelacqua,  and  J.  B.  Schneider,  “Backward- Wave  Meta-Materials  for  Perfect 
Eenses,”  IEEE  AP-S  International  Symposium  and  URSI  Radio  Science  Meeting,  San  Antonio,  TX,  Jun. 
2002. 

[6]  M.  W.  Eeise,  P.  J.  Bevelacqua,  and  J.  B.  Schneider,  “Effects  of  Surface  Waves  on  the  Behavior  of 
Perfect  Eenses,”  Physical  Review  B,  vol.  66,  no.  3,  035 113  (five  pages),  2002. 


7 


